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Cn 1 Abstract. In this note we consider two different singular limits to hyperbolic system of conservation 

f^ / laws, namely the standard backward schemes for non linear semigroups and the semidiscrete scheme. 

fl^ ' Under the assumption that the rarefaction curve of the corresponding hyperbolic system are straight 

rf\ , lines, we prove the stability of the solution and the convergence to the perturbed system to the unique 

solution of the limit system for initial data with small total variation. 

Oh ; S.I.S.S.A. Ref. 85/2000/M 

j^ , 1. Introduction 

Consider a hyperbolic system of conservation laws 

ut + f{u)x = 



(1-1) ^ (n \ f \ 

CN i where u G R" and / is a smooth function from an open set ft C R" with values in M". Let Kq be a 



compact set contained in il, and let 6i sufficiently small such that the compact set 
(1.2) ilTi = |u e R" : dist(M, Kq) < Si} 



is entirely contained in fl. 

We assume that the Jacobian matrix A — Df is uniformly strictly hyperbolic in i^i, i.e. 



o 

0^ 
O 
O 
O 

^. 

C^ . (1.3) mini Xj{u) - \^{u)} >c> 0, yu,v e Ki 

where we denote by A.; the eigenvalues of A, Xi < Xj. Let r^, T be the its right, left eigenvectors. 

In this setting it is well known that if ?io(— oo) G Kq and Tot.Var.(uo) is sufficiently small, there exists 
^^ ! a unique "entropic" solution m : [0, +oo) ^ wo + i^(M,R")nBV(R,M") in the sense of [|| . Moreover these 

$-H ' solutions can be constructed as limits of wave front tracking approximations and they depend Lipschitz 

^ continuously on the initial data. 

For a special class of systems, called in [|6| Straight Line Systems, i.e. systems such that 

(1.4) (Dn)n = 0, 



very recently it has been proved that solutions to (1.1) can be constructed as L^ limits of solutions to 



different singular approximation of the hyperbolic system: 

• Vanishing viscosity approximation |^ . This is the limit as e —* of the solution u' (i) of the system 

ut + f{u)x - e-u^x = 0. 

• Relaxation approximation |S. H. While in case 1) the perturbation is parabolic, in this case we 
consider a hyperbolic perturbation, namely 

"t + f[u)x = e{K^Uxx - utt), 
where A is strictly bigger than all the eigenvalues of Df{u). 
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Godunov scheme 1^. This is a discrete scheme obtained from (LI) by considering differential ratio 
instead of derivatives: 

u{n+ f,j + l) =u{n,j + 1) + ^ f{u{n,j)) - f{u{n,j + 1)) 

where for stabihty condition it is assumed that < Ai < • • • < A„ < Ax/ At. 

The main idea behind these approximations is to obtain uniform BV estimates for t > 0, if the initial 
data Wo are of sufficiently small total variation. 

This task is achieved by decomposing the equations satisfied by Ux, or Ux and uj in |y|, or u{n,j) — 
u{ti,j ~ 1) in gl, as n scalar perturbed conservation laws, coupled by terms of higher order. These 
terms are then considered as the source of total variation. For the special case of straight line systems, 
a decomposition of Ux which makes the source terms integrable is the projection along the eigenvectors 
r, of the Jacobian Df{u): 

Ux = y^fVi. 

i 

Once it is proved that the L^ norm of the component w' is bounded, by Helly's theorem there exists a 
subsequence u'^'^ converging to a weak solution u{t) of ( |1.1| ) as fc — » oo. 

To prove the uniqueness of the limit u{t), one consider the equation for a perturbation h of the singular 
approximations. We observe that h — Ux 'is & particular solution of such system. A generalization of the 
arguments used to prove an a priori bound on the total variation of u shows the boundedness of the L^ 
norm of the components /i', where 

h ~ 'S^ hWi{u). 

i 

By a standard homotopy argument Isl, this yields the stability of all solutions of the approximating 
system. Since the Lipschitz continuous dependence on the initial data is uniform w.r.t. both e and t, in 
the limit we obtain a uniform Lipschitz semigroup. 



Finally it is well known that a uniform Lipschitz semigroup of solutions to (1.1) is uniquely defined if 
we know the jumps conditions of the entropic shocks, see Q. In this case, because of the condition (1.4) 
and because in the scalar case the solution u"^ converges to the entropic solution, an argument similar to 
the one in H implies that the jump conditions coincide with the scalar jumps along the eigenvectors r^. 

Thus, under the assumption ( |l.4| ), the limit semigroup is independent on the approximation and 
coincides with the solution constructed by wave front tracking using the classical Lax Riemann solver. 

Li this note we want to extend the previous approach to the following cases: 

1. semigroup approximation [pi p|. This is obtained as limit of the system 

u(t,x) — u(t — e,x) ,, , ,x , , 

-^■^-!- i '—!- + A{u(t, x))u(t, x)x = 0. 

e 

This is the standard backward scheme for non-linear semigroups. 

2. Semi-discrete schemes |l), i.e. infinite dimensional ODE defined by 



d_ 
di 



u{t,x) + -(^f{u{t,x)) - f{u{t,x~ e))) = 0. 



We will prove that as e — > the limits of the respective solutions converge to a unique solution to (1.1), 
and that this limits defines a Lipschitz continuous semigroup S on the space of function with small TV. 
Moreover this semigroup is perfectly defined by a Riemann solver which, as explained above, coincides 
with the classical one. 

The same can be proved for quasilinear systems as in p, pi , but for simplicity we consider only systems 
on conservation forms. 

Without any loss of generality we assume that 

(1.5) min{Ai(w)} = k > 0, max{Ai(u)} — K < 1, 

for all u in the compact set Ki. The second condition is needed only in case 2). 
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2. Approximation by semigroup theory 



We consider in this section the case 1) of Section [|, i.e. the foUowing singular approximation to system 
of conservation laws: 

u(t,x) — u(t — e.x) ,, , ,^ , 

(2.1) ^^^ f ^ + A(u(t,a;))u(t,a;), = 0, 

where we recall that u G M" and A[u) — Df{u). By the rescaling t -^ t/e, x -^ x/e and setting for 
simplicity Un{x) — u{n, x), we obtain the evolutionary equations 

(2.2) U„ - Un-l + A(un)un^x = 0. 

It is easy to prove that if the BV norm of u„_i is sufficiently small, then u„ exists: in fact the solution 
can be represented as 

Un{x) ^ / exp<^ / A^'^(un{z))dz'> A^^{un{y))un~i{y)dy, 

and since the eigenvalues of A are positive we have that 

||u«||^ < CTot.Var.('u„_i), 
C being a uniform constant of ||A^^||oo in the compact set Kq. 

2.1. Projection on rarefaction curves. We now start the procedure explained in Section nl By 
projecting the derivative along the eigenvectors r.i(w„) of A{un) 

(2-3) Un,x^^vl^ri{Un) ^^vl^n^n, 

i i 

the equations for the components v^ are 

i i i 

This can be rewritten as 

i i i,j 

The left-hand side is in conservation form, and we consider the right-hand side as the source of total 
variation. If we assume as in the introduction that {Dri)ri{u) = 0, the function ri[u) — ri{v) is zero when 
u — V is parallel to ri{u) = ri{v). Thus we have 

(2.5) n{u) -n{v) = y^^aj{u,v){P{u),u-v), 

where aj{u,u) = rj{u). Using (pT^), the expansion (p]j) thus becomes 

= H H'^kinK-A + Y. K]k{n)vivl. 
To estimate the source terms in (|2.6|), we first consider the case of two linear equations. 



2.2. Analysis of the linear case. Consider a single linear equation 

(2.7) Vn - Vn-l + XVn,x =0, A > 0. 

We can find the fundamental solution to the previous equation by means of Fourier transform: we have 

Jr 
and substituting 

c(n,e)-c(r^-l,O-^A^c(n,O-0 =^ c(n, = -^^^rrr- 

(1 - iXQ 
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r^. 



Figure 1. Integration path on the complex plane, where P = i{X — /i)/A/x. 



In particular the fundamental solution has co(C) = l/27r, so that 

(2-8) ^"*'^'' " A V A/ (-^^_;^\| X[o,+oo)(a:)- 



Consider two equations of the form (2.7), 

(2.9) Vn - Wn-1 + XVn,x = 

^n ^n— 1 ~r ^Zn,x — U 

with initial data va{x) — S{x) and zo{x) = d{x — xq), and assume without any loss of generality that 
A > /i > 0. We can compute the intersection integrals: denoting with d{n,^) the Fourier coefficients of 
Zn{x) we have 



(2.10) 



Y, / vr,{x)zn{x)dx = ^2^ / c{n,Od{n,~^)e-'i-°d^ 



1 /■ ^ 1 



-i(xo 



dt 



If ^ is considered as a complex variable, we can let N — > +cxi only in the region where 

Z 



{eeZ:|(l-zAe)(l+*MOI<l} 



i.e. outside the regions depicted in Figure H. Deforming the path to avoid the region Z, we can pass to 
the limit: 



+ C30 



n=0 



y^ / Vn{x)Zn{x)dx 



1 



-i^Xo 



-de 



2^^ A 1 (i-»ac)(i+(pC) 

1 r j^^aoo+M^.-'C-o 



27r / (l-iAC)(l +1^0-1 



de 



By means of complex analysis we have finally that 



+ 00 



(2.11) 



n=0 



Pi^o) = X! / '"nix)Zn{x)dx = 



/{\- fJ.)-CXp({\- ^)/{\^)xQJ Xq < 



.l/(A-/i) 



xo > 



In fact, depending on the sign of xq, the integration along the line 7 is equivalent to the integration 
around the pole or the pole P = i{X — /i)/A/i. 
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2.3. BV estimates. Now to prove that ( ^.2| ) has a solution with uniformly bounded total variation. 
Define the functional 

(2.12) Q{n) = Q{u„,Un-i) ^Yl / Po{x - y)\^\vlix)vi{y)\ + \vl_-i^(x)vi{y)\ + \v!^{x)v'„_j^{y)\jdxdy, 

where P is computed substituting to A — /^ the constant of separation of speeds c, and taking the minimal 
value of the exponent (A — fi)/Xii: 



Pn{x) 



1 



/c-exp\c/{K{K -c))xo] xq 



<0 



.1/c 



a;o >0 



We recall that c and Kappa are defined in the introduction. 
Using the same analysis of [g| , we see immediately that 



(2.13) g(n)-Q(n-l)= ^ / Po{x-y){\vl{x)vi{y)\-\vl_,vi_A 

+ Y, f Po{x-y){K.i{x)vi{y)\-\vl^_^vi_^\} 

+ Z! / ^o(a;-y){K(a;)w^_i(2/)| - \v\_^vl_^\^dxdy 



< — (1 — C max Tot.Var.(Mm) 






where C is c constant depending only on i?,'^,, -^.j^,, k, K and c. Thus if (5o is sufficiently small, using ( 2.13 ) 
we have 



Tot.Var.(wi)+CoQ(ui,uo) < (5i and — |Tot.Var.(u) + Cog(u)| 



<0, 



where the constant Co is big enough, independent on bi^. This proves that the solution m„ has uniformly 
bounded total variation for all n e N. 



2.4. Stability estimates. We now consider the stability estimates of (2^). The equations for a pertur- 
bation M + (5/i as (5 ^ are 

(2.14) hn - hn-l + {A{Un)hn)^ = (L'^(w„) U„^a;) /l - (DA(un) hn)un,x- 

Using the same projection of ( |2.3| ), i.e. 

i 

we have that the equations for the components /ijj are 

(2.15) hi, - /ij,_i + (A^,„/ij,)^ = J2{^k,nK_iV^{ll,,aj{u„,Un-i)) - hl,vi{ll^, (L>r,,„)rj,„)j 

= J2 H{n)K-i< + E ^'(")'^>«- 



Using the same analysis of , it is easy to prove that a functional as in Section p.3| gives the stability of 
the solution. 
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3. Approximation by semi-discrete scheme 

We now consider the case 2) of Section |], i.e. the foUowing singular approximation to system of 
conservation laws: 

(3.1) ^w(t, x) + i (/(w(t, x)) - f{u{t, x-e)))^0, 

where u G M". By the rescaling t — > t/e, x -^ x/e, we obtain the evolutionary equations 

(3.2) Unit) + f{un{t)) - /(Un-l(t)) = 0. 

The equation for the "derivative" Wn = u„ — m„_i are 

(3.3) Vn{t) + /K(i)) - 2/(u„_i(i)) + f{un^2{t)) - 0. 

3.1. Projection on rarefaction curves. The vector Vn is now decomposed along the eigenvectors ri^n 
of the Riemann problem u„_i, u„: we have 

Un{t) + ^ XijiVl^Ti^n = 0, 
i 

z 
i,3 i,j 

where Ai_„ and r^^^ are the eigenvalues and right eigenvectors of the average matrix 

A{un,Un-l) = j Dfiun-l + (u„ -U„-l)sjds. 

If we assume the condition ( |1.4D , the functions {Dri^n)fj,n and r.i^n — i'i,n-i are zero when u„ — m„_i 
and Un-i — Un-2 are parallel to ri_„ = ?'i,n-i- Thus we have 

(3.5) rj,n»n,n= '^aj^nvi, 



as in Section 2.1. Using (p.5|), the expansion (pA) thus becomes 



(3.6) < + A,,„<-A,,„_i<_i =Xi/„(t)l;>,'; + XG„(i)^'>Ll• 



To estimate the source terms in (3.6), we consider the case of two linear equations. 



3.2. Analysis of the linear case. Consider a single linear equation 

(3.7) < + A< - A<_i =0, A > 0. 

We can find the fundamental solution to the previous equation by means of Fourier transform: defining 
the periodic function 



c(t,x) = X«„(i)e™^ 
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Figure 2. Integration path on the complex plane, where P = fi/X. 



we have that the equation satisfied by c is 



ct = ^^„e™^ = A^(i;„_i - i;")e*"^ 

n n 

= A(e" - l)c, 
whose general solution is 

c{t, x) = c(0, a;)exp(A(e*=^ - l)t] . 
In particular the fundamental solution starting at hq has c(0,x) = expiiuQx), so that if tiq — 

'O X <0 



(3.8) 



Vnit) = 



{Xty'/n\-exp{~Xt) x>0 
If now we consider two equations of the form (pj), 
(3.9) Vn + A(v„ - W„-l) = 0, 

Z„ + H{zn - Zn-l) = 0. 

we can compute the intersection integrals: denoting with d{t,x) the Fourier transform of Zn{t) and 
assuming that A > /i > 0, we have 



(3.10) 



n— — oc 



/ V Vn{t)zn{t)dt ^ / — / c{t,x)d{t,-x)e~™<'''dxdt 

Jo „__ Jo 27r7o 

)fA(e"-l)i + Ai(e-"-l)n< 



exp 



^0 

27T 



1 

2^ 

1 
27r 7o A(e" - l) + /i(e"" - l) 

1 f z-"" 

2^7^ (z-l)(Az-/z) ^' 



"dxdt 



dx 



where 7 is the path represented in Figure y. 
By means of complex analysis we have that 



(3.11) 



Pino) 



/ X! V.a{t)Zn{t)dx = < 

•^0 n=-oo V 



1/(A - /i) . (A/Ai)"° no<0 
/(A - /i) no > 
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3.3. BV estimates. Now to prove that (3J) has a solution with uniformly bounded total variation. By 
defining the functional 

+00 
(3.12) Q{u{t))^J2 E P{n-m)[vi,{t)vUt)+vl,__,{t)vUt)+vlMyLMt)}dxdy, 

i<Cj n,m— — oc 

where P is computed using the constant of separation of speeds c instead of A — /i and 1 + c/ K instead 
of A//i, since the left hand side of ( |3.6D is in conservation form, we conclude immediately that 

Tot.Var.(u(0)) + CqQ{u{o)) < Si, — |Tot.Var.(M) + CoQ{u)\ < 0. 

This concludes the proof of bounded total variation. 



3.4. Stability estimates. Finally we consider the stability estimates of (3.1). The equations for a 
perturbation u + <5/i as (5 — > are 

(3.13) hn{t) + Df{un)K - Df{un-l) = 0. 

Considering the projection 

i 

we have that the equations for the components /i^ are 

At this point it is clear that a functional as in Section p.3| proves the stability of the solution. This 
concludes the proof. 
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